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CO ' By means of a cosmographical approach to the redshift drift, together with a series development 

of the Hubble parameter for an arbitrary f(R) theory, we obtain a restriction on / and some of its 
derivatives evaluated at the present epoch. It is show that, when applied to specific forms of f(R), 
this restriction (which depends on the values of the kinematical parameters today) leads to bounds 
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on the parameters of the theory that are more stringent than those previously obtained. 
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I. INTRODUCTION 
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The interpretation of several sets of data (such as those obtained from type la supernovae, large scale structure, 
5^ ' baryon acoustic oscillations, and the cosmic microwave background) in the framework of the Standard Cosmological 
■ Model (SCM) (based on General Relativity (GR) and the Cosmological Principle (CP)) indicates that the universe 
is currently undergoing a phase of accelerated expansion. The most commonly accepted candidates to source such 
an expansion are the cosmological constant and some unknown type of matter dubbed "dark energy" (for a complete 
list see [l|). While the energy density associated with the cosmological constant that is inferred from astronomical 
observations is approximately 120 orders of magnitude lower than the value predicted by field theory (see for instance 
Q), the scalar field used to model dark energy has features that are alien to those displayed by the scalar fields of 
\q \ particle physics [jj. An alternative to describe the accelerated expansion is to assume that a theory different from 
GR governs the dynamics of the gravitational field for very low values of the curvature. The simplest generalization 
f^*) . of the Einstein-Hilbert Lagrangian is given by 
CO 



s = J V=gf(R)d 4 x, 



where the dependence of the function / on the scalar curvature R is to be determined by several criteria (such as 
matter stability Q| , absence of ghost modes in the cosmological perturbations Q , and correct succession of cosmological 
eras Q), and by the comparison of the predictions of f(R) theories with observation 1 (for instance, the stability of 
cosmological perturbations Q). Several forms for f(R) have been constructed in order to successfully satisfy these 



constraints (for instance those given in jld lllj). 

Currently a lot of effort is devoted to the study of quantities and effects that have the potential of discriminating 
between different f(R) models, and between these models and the SCM. Among them, we can mention the growth rate 



of matter density perturbations (see for instance 12[), the enhanced brightness of dwarf galaxies |13| . the modifications 



of the 21cm power spectrum at reionisation 14(, the specific angular momentum of galactic halos [l5j, and the number 
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1 For reviews about different aspects of f(R) theories, see @, 01, SB- 
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counts of peaks in weak lensing maps [3] . We shall show that another observable that can be of use along these lines 
is the redshift drift (Az/At), that is, the time variation of the cosmological redshift caused by the expansion of the 
Universe. The redshift drift was first considered by Sandage [l7j], and the effects of a nonzero cosmological constant 
on it were presented in 18j. In 1998, Loeb 19( discussed the possibility of actually measuring the redshift drift 



of cosmological sources using spectroscopic techniques available at that time. Since then, the potential of planned 
extremely large observatories, such as the European Extremely Large Telescope (E-ELT), the Thirty Meter Telescope 
and the Giant Magellan Telescope, and the development of ultra-stable spectrographs with high-resolution such as 
the COsmic Dynamics and EXo-earth experiment (CODEX) proposed for the E-ELT, toghether with simulations, 
showed that the measurement of the redshift drift is feasible in the near future (for a review, see ) , leading to what 
has been dubbed "real-time cosmology" in 
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When compared to other cosmological observables, the redshift drift has the advantage that it directly tests the 
dependence of the Hubble parameter with the redshift, hence probing the dynamics of the scale factor. Another 
feature of this observable is that it does not depend on details of the source (such as the absolute luminosity), or 



on the definition of a standard ruler. As demonstrated in |22|, the measurement of the redshift drift would allow 



the test of the Copernican Principle, thus checking for any degree of radial inhomogencity. This issue was further 



discussed in [231 ]. where it was shown that the redshift drift is positive for sources with z < 2 in the ACDM, while in 
Lemaitre-Tolman-Bondi models (see for instance 24 1 ) is negative for sources observed from the symmetry center [2o] | . 
The redshift drift can also be used to constrain phcnomcnological paramctrizations of dynamical dark energy models 
(see 0-01). In this work, we shall extend the applications of the redshift drift by showing that it can be employed 
to set limits on the parameters of an arbitrary f(R) theory. This goal will be achieved by combining a cosmographical 
approach to the redshift drift with the corresponding dynamics. 

Cosmography is a mathematical framework for the description of the universe, based entirely on the CP, and on 
those parts of GR that follow directly from the Principle of Equivalence [3]. It is inherently kinematic, in the sense 
that it is independent of the dynamics obeyed by the scale factor a(t). The cosmographical approach has already led 
to interesting results in the framework of f(R) theories [IB-El- It will be used in Sect. [TTJ to obtain a series expansion 
of the redshift drift in terms of the cosmological redshift z and the derivatives of the scale factor. Since the redshift 
drift depends on the explicit form of the Hubble parameter H(z), we shall show in Sect. |Hl] that its series expansion 
can actually be written in terms of f(R) and its derivatives, using the Equations of Motion (EoM) for the f(R) theory 
in the metric version. By comparing the two series, we show that there exists relations between f(R), its derivatives, 
and the kinematical parameters, which impose constraints on the parameters of a given f(R). By using one of these 

and the exponential gravity theory introduced by 



relations in two examples: that proposed by Hu and Sawicki 
Lindcr 
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34j. In both cases, we find new limits for the parameter space of each theory and compare our results with 



previous constraints. We close in Sect. [TV] with some remarks. 



II. A COSMOGRAPHICAL APPROACH TO THE REDSHIFT DRIFT 



In this section we shall sketch the calculation that leads to the expression for the series expansion of the redshift 
drift in terms of z in homogeneous and isotropic space-times. The redshift of a photon emitted by a source at t = t s 
that reaches the observer at t = t a b s is given by 



Z(t s ) 



a(t b 



a(t s ) 



1. 



(1) 



The time variation of the redshift is obtained by comparing this expression with the one corresponding to a photon 
emitted at t = t s + At s , that is Az = z(t s + At s ) — z(t s ). To first order in At b s and At s , it follows that [19 | 



Az 

At , 



a(t obs ) - a(t s ) 
a(t s ) 



(2) 
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Using the definition of z and H(z) in this equation we get the expression of the redshift drift in terms of z, namely 

Az 



M , 



= (l + z)H obs -H(z). 



(3) 



In the case of the SCM, the redshift drift can be written as a function of the cosmological parameters using the 
corresponding expression for H(z) as follows fl9| : 



(Az/ At) = H (l + z)- (O ro , (l + zf + n r , (l + z) 4 + f2 A)0 ) 



1/2 



(4) 



where the subindex means that the corresponding quantity is evaluated at t = today. As discussed in the Introduc- 
tion, we would like to pursue here a mixed approach partially based in the concept of cosmography (see for instance 
35|). Cosmography consists in constraining the history of the scale factor (and a posteriori, of quantities defined 
in terms of a(t)) from observational data, without using the dynamics obeyed by the gravitational field. The series 
development of the scale factor around to is given by 



a(t) = a 



1 + H (t - t Q ) - ^q H 2 (t - to) 2 + ijofloC* - M 3 + ^s H 4 (t - to) 4 + 0([t - t } 5 ) 



(5) 



where the so-called kincmatical parameters arc defined by 



„, . 1 da , . 1 d 2 a 



1 da 

a(t) dt 



1 d 3 a 



1 da 
i(t) dt 



1 d 4 a 



1 da 
a(t) dT 



In order to use Eq. ([5]) for the calculation of the redshift drift, we need to express t s in terms of known quantities. 
This can be done through the physical distance travelled by a photon emitted at t s and observed at to, given by 



D = cjdt = c(t - t s ). 

This equation, toghether with (TT|) yields 

H D , 2 + qoH^D 2 6(1 + q a ) + jo H^D 3 



z(D) = 



which can be inverted to 



D(z) 



Ho 



i + - )z+ n +qa + 2 



% Jo 



G 



1+ 2?o(l + go)+ g?o 



O 



HoD 



-I 4> 



5 . s \ 3 4 . 



(6) 



(7) 



(8) 



Using Eqs. (|6]) and we can write a(t s ) as a power series in z, with coefficients that are functions of the kincmatical 
parameters. After some algebra, it follows that 



= -Hoqoz + ^H Q (ql - jo) z 2 + ^H 



2 (so + Mojo) + jo -ql-ql 



z 3 + 0(z 4 ). 



(9) 



This equation gives the expression of the redshift drift up to the third order in the redshift of the source, in terms of the 
value of the kinematical parameters at the present epoch. Note that the coefficient of z n involves the (n + l)-derivative 
of the scale factor. 

In the next section we shall obtain an expression for the redshift drift using the dynamics of an arbitrary f(R) 
theory toghether with the development in terms of z. 

III. THE REDSHIFT DRIFT IN f(R) THEORIES 



Let us recall that the redshift drift can be expressed in terms of H(z) as follows 

Az 



At 



(1 + z)H - H(z) 



(10) 
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The dynamics of a given theory establishes the form of the function H(z). In the case of the SCM, we obtain for the 
redshift drift the expression given in Eq. ^ using Friedmann equation. In the case of f(R) theories, the EoM that 
follow from the variation of the action 



S = J dW^[/OR)+A„atter] 



w.r.t. the metric must be used instead. For the FLRW metric and considering a prcssureless cosmological fluid, these 
equations are (see for instance, [3]) 

f'R - 2/ - 3/" + - 3f'"R 2 = T, (11) 

f'Ru + \f + V"^=T tu (12) 

• 2 

where dot and prime denote, respectively, derivative w.r.t. t and R, Ru = 3'd/a, R = 6(^ + %), and T is the trace of 
the energy-momentum tensor. From these, an equation for H can be obtained, and it reads 

H = --^p {6H 2 f - 2 Pm - f + Rf>) . (13) 

Replacing Eq. (fTO)) , we obtain 



A, = do^) = lf(R)-Rf(R) + ^ 2 f { R)-2p m 

At y ' a{z) 1 ; 6 R(t)f"{R) 

where R, H, and p m are functions of z. Expanding this equation in powers of z and using the definitions of the 
kinematical parameters, the redshift drift up to first order in z for an arbitrary f(R) is given by 

^(z) = { [ Jo 2 (l - q ) + 4(1 - j ) + 2q OJO (l + go) + g 2 ('3 - go)] 36iJ 4 (/^') 2 

+ [(so + ql + 6 + 8<fo)(/ - 2p m ) + (2j (l + go) + qo(s + <Zo + 6%) - 4) 6ff 2 / ] ft 
+ [Qg Hlf - 2 Pm + f ] [(q - Jof + 4(1 + g - Jv)] GH 2 /'"} 7^773 + 0(z 2 ). (15) 

where we have used the expressions R = 6(H + 2H 2 ), H = —H 2 (l + g), and H = H 3 (j + 3q + 2). The second order 
term involves derivatives of the scale factor higher than the fourth. 

Comparing expression (fT5|) with the linear term of Eq. ([9]) we find a relation between f(R), its derivatives and the 
kinematical parameters, all evaluated at t =today: 

[(so + ql + 6 + 8gb)(/ - 2p m ) + (-4 + 6q 2 + 2j a + q s + 2j q + ^)6i/ 2 / ] ft+ 

+ [(9o - Jo) 2 + 4(1 + qo - jo)] [(/o - 2p m + 6 qo H 2 .Qft' + 6H 2 (ft) 2 ] 6H 2 = 0. (16) 

Since we are working at first order in z, this equation is a necessary condition for any f(R) theory to describe the 
variation of the redshift drift with z. By equating higher orders of z from Eqs. ([9]) and (|15j) we would obtain more 
(actually, an infinite number of) necessary conditions on f(R) and its derivatives. If the theory under discussion is to 
describe the redshift drift at all orders in z, all these conditions should be satisfied. Let us remark that the condition 
that would be obtained from the second order in z would not be restrictive, since it will depend on the value of the 
kinematical parameter £ (proportional to the fifth derivative of the scale factor) , for which there are no observational 
limits available. 

We shall see next how Eq. (|16p constrains the value of the parameters of a given f(R), by applying it to two 
particular forms of f(R). 
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A. Example 1: The Theory of Hu and Sawicki 



Let us start with the theory introduced by Hu and Sawicki [llj , which is given by 



f(R) = R-m C2{R/m2)n + 1 , (17) 

where n > 0, c\ and C2 are dimensionless parameters and the mass scale is m? = k 2 pq/3, with po the average density 
today. 

For values of the curvature high compared with m 2 , f(R) may be expanded as 

lim f(R) ~ R - —m 2 + %m 2 (m 2 /R) n (18) 

m 2 /_R->0 C2 Cn 



and, at finite Ci/c 2 ., the theory can approximate the expansion history of the ACDM model 11] . In this regime, the 
parameters c\ and C2 must satisfy the relation 

— «6^. (19) 
For the flat ACDM expansion history, Eq. (fTg|) yields 

i? ~ 3m 2 ( a - 3 +4fi^) , (20) 

i?o - m 2 [^--9] , (22) 

/o - 1 - «| f jp- - 9) • (23) 

Using Eq. (|23|) . ci/c 2 , can be expressscd in terms of n, /q and tt m .o- In addition, higher order derivatives of / can 
also be written in terms of the same quantities. Hence from now on we set f2 m> o = 0.274 ± 0.007 and leave /q 
and n the only free parameters of the theory. With these considerations, Eq. (|16|) yields 



and at the present epoch, 



/o(n) = j { [q {8 + 5q - 6j + (go - jo) 2 ) + (jo - 2) 2 ] An + [70(70 - 2g - 4) - 14 - 2g o (10 + go) - 3s ] 6tt m , 

+ [93(8 + go) + 6 + so] 9f) 2 n , + 4 [2 (g (10 - 2j + 3g ) + 6 + s ) + g (go - jo) 2 ] } , (24) 

with 

A = [q (8 + 5g - 6j + (go - Jo) 2 ) + (jo - 2) 2 ] An - [2gt>(13 + 2g ) - jo(jo - 2) + 18 + 3s ] 6fi m>0 

+ [g (8 + g ) + 6 + s ] 9f) 2 m , + 4 [2(10 + 5g 2 + 16g - 2j a (2q a + 1) + s ) + g (g - jo) 2 ] . 

Expression (|24[) gives a relation between the parameters /g and n in terms of fi rn .o and the kinematical parameters 
H, q, j, and s, all evaluated today. We shall take the values H$ = 74.2 ± 3.6 km/s/Mpc 38[, go = —0.669 ± 0.052, 



jo = 0.284±0.151 and sq = — 0.680±0.456 [3l|. In Figure[T]we plot the relation f' (n) provided by our cosmographical 
approach to the redshift drift combined with the expansion of the expression of H(z) for f(R) theories. The curve 
tends asymptotically to /q = 1, which corresponds to the GR limit. The plot also displays the limit obtained from 



solar system tests, given by |/q — 1| < 0.1 [ll|. We find that actually 1 — /q < 0.1 and, from this limit, values 
for n larger than 12.8 are favoured. This result improves the previous bound for the parameter n obtained from 
measurements of the linear growth of structure and SNIa luminosity . 
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redshift drift constraint 
|f„-1|<0.1 
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FIG. 1: Constraints on the parameter space of the theory by Hu and Sawicki. The red line corresponds to the 
relation between /q and n that follows from Eq. (|16p . The dashed lines represent the error propagation arising from 
the kinematical parameters and the density matter at the present epoch. The shady region indicates the values of /q 

in agreement with solar system tests. 



B. Example 2: Exponential Gravity 

Next we shall analyze the restrictions that follow from Eq. (|16|) on the choice of f(R) proposed by Linder (j^ : 

f(R) = R-cr(l-exp(-R/r)), (25) 

where c and r are two (positive) parameters of the model. This f(R) was specifically designed to (i) avoid the inclusion 
of an implicit cosmological constant (since it vanishes in the low curvature limit), (ii) reduce to GR for high values of 
the curvature, (iii) incorporate a transition scale (given by r) to be fitted from observations (instead of set equal to 
i?o), and (iv) restore GR for locally high curvature systems such as the solar system or galaxies. As shown in 34j . 
the product cr is given in terms of f2 m by 

cr = 6m 2 (Q" 1 - 1). (26) 

The use of Eq. (|16[) for the current choice of f(R) yields a relation between the dimcnsionless parameters c and r/m 
given by 

c (^2) = ~b { _6 ^2^ m ^ q °^ + *?o) + 6 + s ) fi m - (2(1 +j qo + qo+jo) + s - $)] 



+36 [A(q + 1 - j ) + (q a - j ) 2 ] (n m + 1)} cxp ( 6(1 qo) 



(27) 



with 

B = I [q {8 + qo) + 6 + s ] 



cxp 



/6(l-g ) 



+6 |g (9o + 7g + 4 + so) + 3o{jo - 2) - [4{q + 1 - j ) + {q - j ) 2 ] cxp (^L^l^j | 



-36 [4(1 + 2q Q ) + q 2 (5 + q ) + j 2 Q {l + q ) - h(2q 2 + 6q + 4)] , 

which is plotted along with Eq. (|26j) in Figure [21 The strips formed by both curves and the corresponding errors 
juxtapose for all values of c > 12, which imply that r/m 2 < 1.4. Notice that the range of possible values for the 
parameter c that follows from our method considerably improves the previous bound (c > 1.27) obtained in [4o| . 
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FIG. 2: Constraints coming from our approach to the redshift drift (continuous line) and the relation ([26]) between 
the parameters of the theory (dashed line). In both cases, the thin lines represent the error propagation arising from 
the kincmatical parameters and the matter density at the present epoch. 



We have developed a cosmographical approach to the redshift drift and compared it with a series expansion in 
terms of z of the Hubble parameter for an arbitrary f(R) theory. As a result, a relation between f(R), its derivatives 
and the value of the kinematical parameters, all evaluated at the present epoch, was obtained. This relation must 
be satisfied for a given f(R) theory to describe the redshift drift at first order in z. The comparison of higher-order 
terms would yield relations between the quantities mentioned above, that involve kinematical parameters of order 
higher than the fourth (which have not been determined with enough precision yet). It is worthile noting that this 
method does not depend on the actual measurement of the redshift drift which, as mentioned in the Introduction, is 
not available yet. 

We have shown that the relation given by Eq. (fT6j) can be used to set limits on the parameter space of a given f(R) 
theory. As examples, we have selected the theory proposed by Hu and Sawicki, and the exponential gravity theory 
introduced by Linder. For the former, we found a constraint on the possible values of the parameter n that improves 
previous bound obtained from other cosmological data. In the case of the latter theory, our method ameliorated the 
previous bound for one of the parameters of the theory (c), and provided the first bounds on the second (r/m 2 ). 

It is worth noting that the approach presented here, which is based on the redshift drift in f(R) theories, can in 
principle be applied to other observables quantities relevant to cosmology, as well as to other alternative theories of 
gravity. 

To close, we would like to emphasize that the bounds obtained by the method developed here can be analyzed 
toghether with those coming from the observations mentioned in the Introduction as well as other means (such as 
energy conditions |41(), with the aim of deciding whether a given f(R) theory is consistent with the available data. 



IV. DISCUSSION 
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